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1. Introduction 

Following the pioneering papers of Ref. [|], it has been realized J2], |3| that large rapidity gaps (LRG) diffrac¬ 
tive processes are suppressed due to the rescattering of the spectator partons. The calculations of the 
resulting gap survival probability (SP)@, |6], [|, ||, ||, |j, [nj have since then been at the focus of high energy 
phenomenology. The incentive for this activity has been the need to obtain reliable SP estimates for Higgs 
diffractive production at the LHC. The analogous process of diffractive hard di-jet production has been 
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Figure 1: The contaminating soft rescatterings between partons, that produce secondaries which may fill the rapidity 
gaps (Fig. \^-a), and the lego-plot (Fig. [^- b) of our process. The probability that there are no additional 
soft interactions shown in Fig. |J-a, gives the survival probability for the diffractive di-jet production. 


measured at Fermilab ||TT|| and HERA ||l2| and this may serve as a laboratory to check the validity and 
reliability of the proposed calculations. 


SP is the probability to have a simple diffractive colorless final state configuration, in which the LRG 
are preserved, regardless of the strong interaction of the rescattered soft partons (see Fig. [I]). Rs calculation 
depends crucially on soft scattering physics, for which a theory is still lacking. The main phenomenological 
tool utilized in these calculations is the Regge soft Pomeron. In spite of the progress made in understanding 
the soft Pomeron structure within the framework of npQCDfl3|, |l4| |l5| , [T6] , o, m n, n p , we are still 
unable to predict the parameters associated with the phenomenological Pomeronp^l, and to create a theory 
for soft scattering. Consequently, no theoretical approach exists which allows us to calculate a SP value to 
a specific given accuracy. As a substitute, the models which were developed utilize relatively simple soft 
Pomeron parametrizations and the rescattering process is approximated by eikonal-type models. These 
models satisfy the general principles of unitarity, including the Froissart and Pumplin bounds, and allows 
us to analyze the experimental data in accordance with these general principles. 


At present, we can obtain only phenomenological estimates for SP. Although our calculations are 
applicable to any small size colorless system e.g. the production of Higgs mesons, in this paper we refer 
only to the particular case of a di-jet system, as there is experimental data from the Tevatron and HERA. 
In the simplest approach we have three steps in calculating the central hard exclusive production of a di-jet, 
N + N N + LRG + J J + LRG + N, which is illustrated in Fig. [I]. 

1) A description of the rescattering soft interactions. In the eikonal model, the elastic high energy 
amplitude is pure imaginary at high energy, and can be written in impact parameter b-space in the form 
(see, for example, Refs.||, |(|) 


aei(s-h) = z (1 - e-^M)), 

(1.1) 

G in (s;b) = 1 - 

(1.2) 
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s = W 2 , where W is the c.m. energy of the initial hadronic collision, and b is its impact parameter. For 
the opacity f2, we use the exchange of a soft Pomeron with a t = 0 linear trajectory intercept of (1 + A) 
and a slope a’ p . Accordingly, 

/ \ A 


fl(s, b) 


^{tJ S(s ’ b) - 


(1.3) 


S(s,b ) denotes the b-space normalized soft profile function satisfying 



(1.4) 


In the framework of the eikonal model <to, So and the form of the profile function are phenomenological 
inputs, extracted from the experimental data. 

2) We need to find the impact parameter behaviour of the two hard Pomerons, that produce the di-jet in 
the central rapidity region, as is shown in Fig. [fj-a. We denote this observable by . We can use the 
HERA DIS data to find this dependence. In this paper we will use the approach developed in Ref.[^], in 
which the b-dependence of the hard Pomeron has been extracted from the diffractive production of J/’F 
meson in the HERA DIS experiments. 

3) The third step is an actual calculation of the survival probability [§, |b| 


f d 2 b ( da H /d 2 b ) e ~^ s ^ 
f d 2 b ( dau/cPb ) 


(1.5) 


The interpretation of Eq. ( |1.5| ) derives from Eq. (|h^), where the factor e secures that no inelastic 

interactions occur which could change the LRG lego-plot of Fig. 0-b. However, the simple form of Eq. O) 
is valid only for a single channel eikonal model in which only elastic rescatterings are considered. It needs 
to be modified in a more elaborate approach, such as two or three channel models, in which both elastic 
and diffractive rescatterings are included @. 

The above appealing sketch of the eikonal approach demonstrates, as well, the deficiencies which 
are inherent in all presently available models, even those which are more complicated than the eikonal 
rescattering approximation. As it stands, the first and the third steps are stricly phenomenological, and 
only the second step is based on well established pQCD. This shaky theoretical situation has been the 
main reason why we were not active in this field over the past five years. The strong current demand to 
evaluate the LHC diffractive cross sections and their SP, encouraged us to return to this subject. 

This paper has two main goals. To begin with, we wish to calculate a range of possible values for 
exclusive LHC central di-jet (or Higgs) production. To this end, we employ the eikonal model, being 
simple and transparent, so that all uncertainties can be recognized and discussed. We wish to avoid 
more complicated calculations in which it is difficult to separate uncertainties in the values of the input 
parameters, from the qualitative behaviour of the amplitude (mostly as a function of the impact parameter 
b). We shall present and compare results stemmimg from a single channel eikonal model, utilizing two 
profile functions. A constituent quark model (CQM) recently suggested ^ ], and a two channel eikonal 
modeljij. 

Our second main goal is to assess the t-dependence of the calculated SP. The experimental study 
of this differential observable requires, obviously, a significantly higher statistics than presently available. 
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However, as we shall see, the discovery potential associated with this information justifies both experimental 
and theoretical efforts. 


In our approach we use three principles which alleviate our poor knowledge on soft scattering theory. 

1) We aim to successfully describe the soft scattering data associated with our investigation. This includes 
the data base with which we adjust the input information for the construction of the eikonal opacity (or 
opacities) and the model predictions that can be tested experimentally. 

2) Utilize all possible theoretical approaches to control our phenomenology. In particular, the pQCD 
calculation for the production of hard di-jets at short distances. 

3) Even though the LHC is not yet running, we should present our overall LHC predictions so as to be a 
subject for an experimental test even at the preliminary runs, well before the survival probabilities can be 
experimentally assessed. 


For a recent review on SP, see Ref. [241. 


2. Formulation 


In the following we present the main formulae pertinent to our investigations. 


2.1 Basic formulae 


Our elastic t-channel scattering amplitude is normalized so that 


da . „ . . l9 

lu =7T 


( 2 . 1 ) 

( 2 . 2 ) 


Cftot = 4"7r Im f e i (s, 0). 

The amplitude in impact parameter space is given by 

a ei{s,b) = J d 2 qe~ iq ' b f e i(s,t), 

where, t = — q 2 . In this representation we have 

a t ot = 2 J d 2 b lma e i(s,b), 

K,M)I 2 - 

s-channel unitarity implies that | a(s, b) |< 1. When this is written in a diagonalized form we have 

2 Im a el (s,b) = \ a ei (s,b ) | 2 +G in (s,b). (2.6) 

The corresponding inelastic cross section is written 

c Tin = J d 2 b G in {s,b). (2.7) 

s- channel unitarity is most easily enforced in the eikonal approach, where Eq. (jl . 1|) , Eq. (1.2) and Eq. (|1.3|) 


(2.3) 

(2.4) 

(2.5) 


satisfy the unitarity constraint of Eq. (2.6). 
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Figure 2: Central production of two hard jets separated by two large rapidity gaps from the accompanied final state 
nucleons and/or diffracted excited states (both are denoted by N*). Possible rescattering of the initial two 
nucleons (N) is included. 

2.2 The general expression for the survival probability 

The amplitude for central hard di-jet production is shown in Fig. |j. Ah denotes the exchange of a hard 
Pomeron responsible for the production of the two gluonic jets (or Higgs), at short distances. The amplitude 
As includes all possible initial state interactions due to the exchange and interaction of soft Pomerons, 
including the possibility that the two initial nucleons do not interact. 

The expression for the exclusive di-jet production (Fig. |2|) has the form 

A (N + N -> N + LRG + JJ + LRG + N) = (2.8) 

dhard J d 2 k t T h [(p lt - k t ) 2 -,xij T h [ip 2t - kt) 2 -,X 2 j A S (k 2 ) = 

= Ah j <Pb 1 <Pb 2 e i fr t - i '+ i ^AH(b 1 )A H (b2) ((&i + 6 2 ) 2 ) • 

ahard is the amplitude for two hard Pomerons fusion into two jets. This process occurs at short distances, 
and we assume it to be independent of the impact parameters. Ah is the amplitude for the hard Pomeron 
exchange, which is well known in pQCDp5|, [2^ ]. 

The full hard amplitude has a more complicated dependence and can be written as 

Th{Pu x ) = A H {p 2 t)T H (x) —» A H (b)T H {x), (2.9) 
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where the amplitude Th(x) absorbs the energy and transverse momenta dependence of the produced di-jet. 
We do not need to know the explicit dependence of these amplitude on energy and transverse momenta 
for our calculation of SP. 

For the completeness of our presentation we introduce the factor Ah ar d which has been calculated in 
Ref. [27]. This factor is given by 

A hard (HP + HP^JJ) = 


32vr 

~9 


J a s (Q 2 ) ^ r (j){x 1 ,Q 2 )(l)(x 2 Q 2 ) e s ( q2 ’ P t). (2.10) 


HP denotes the hard Pomeron, cf> is the unintegrated structure function and Pp = +P 2 ,i)> where 

Pi t are the transverse momenta of the produced gluons. The Sudakov factor e~ s secures that there is no 
radiation of gluons, with sufficiently large transverse momenta, from the exchanged gluons in Fig. |j. This 
suppression might be cosidered as a contribution to SP stemming from the hard sector. However, we note 
that it is included in the pQCD calculation of the hard process. See Ref. [p?|l for detailed calculations of 
S,(f> and for the kinematical constraints for the process of di-jet production. 

The most important characteristic of the hard Pomeron exchange for the calculation of the SP, is its 
dependence on the impact parameters. This dependence stems from the ^-dependence of the proton-hard 
Pomeron vertex which can be defined in either BFKL||25|1 or DGLAP[26| approach. As we have mentioned, 
this dependence was extracted, in our calculation, from the DIS data^3|, and we shall discuss it below. 

Ag is the amplitude accounting for the initial soft interaction, about which little is known. As stated, 
we start with a calculation of this amplitude in a single channel eikonal model. Accordingly, we consider 
only the elastic rescattering of the two initial nucleons, neglecting the summation over the possible excited 
nucleon states in Fig. |j. 

For this cross section we have 

da 


d 2 p u d 2 p 2 t,dyi dy 2 


(N + N -► N + LRG +JJ + LRG + N) = 


( 2 . 11 ) 


&hard 


J d 2 b 1 d 2 b 2 e i ^^ + i ^^A H (b 1 )A H (b 2 ) A s (fa i + b 2 )‘ 


where pu and yi are the transverse momenta and rapidities of the produced gluons (jets). 

Integrating Eq. ( p.ll ) over pu and p 2 t, we obtain the formula for SP (Eq. ( |1.5| )) in the form 

,1 2 


(I s 


f d 2 b 1 d 2 b 2 A H (bi) A H (b 2 ) A s f (&i + b 2 f 


f d 2 b 1 d 2 b 2 [A H (bi) A H (b 2 )Y 


Let us introduce 


and 


dan 
d 2 b 


= J d 2 b' A 2 H (b-b'^j A 2 H (b'), 


e -n(s,b) _ A 2^ b y 


( 2 . 12 ) 


(2.13) 


(2.14) 


We find that Eq. (2.12) coincides with Eq. CD- We shall show below that Eq. ( |2.14| ) is a consequence of 
the eikonal model. 
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S s (b) = 


p - 


Figure 3: As(b) in the eikonal model as a multi Pomeron exchange chain. P denotes the exchange of a soft Pomeron. 




Figure 4: t-dependence of J/ x V production in DIS. 


In the eikonal approximation the rescattering of two nucleons can be described by the set of diagrams 
presented in Fig. [| To understand the minus sign in front of the first rescattering term, one needs to 


substitute this diagram into Eq. ( 2.11 ), and remember that the single Pomeron exchange, as well as the 
diffractive production of jets, are pure imaginary. This leads to a negative contribution. See, for example, 
Ref.[28]. Summing all diagrams in Fig. [3], we obtain the eikonal formula for 


A s (b) = e -^M), 


(2.15) 


which leads to Eq. ( |2.14 ). It is important to realize, in this context, that different LRG di-jet configurations, 
i.e. JGJ, GJJ and GJJG (which is discussed in this paper), have different A s (b ) and intermediate states 
(see Fig. [| and Eq. (|2.51|) below). Consequently, they have different SP values. 


2.3 A H (b ) 

To fix the impact parameter dependence of the hard amplitude, we return to the model of Ref. |p3|l . The 
t-dependence of the cross section for the diffractive production of J/T is shown in the diagram of Fig. ||. 
The vertex for interaction of the hard Pomeron with the virtual photon leading to the production of J/T, 
can be calculated in pQCD. From this calculation we can find the impact parameter dependence of the 
proton-hard Pomeron vertex. This vertex can be parameterized as an exponential in t 

A H (t) = e -l BGG W, 
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with Bgg = 3.6 GeV |23|]. The corresponding b-space transform of Eq. (p,16|) is 

f / ^ 2 


A H (b) = 


* R 2 h 


exp 


K 


(2.17) 


with R 2 H = 7.2 GeV 2 . In all our calculations below we take this form of An(b). 
Substituting Eq. ( 2.17| ) in Eq. ( 2.13 ) we obtain 


da H 

d?b 


*R 2 h 


exp 


R 2 h 


and 


d 2 b' A h (ft-b') A H {b') = 




_6 2 

R 


H 


(2.18) 


(2.19) 


The uncertainties in the determination of Rh are mostly related to the influence of the excited nucleons 
in the final state (see Fig. j|). To understand better the possible influence of these states on the value of SP, 
we consider the process of diffractive dissociation of the virtual photon to a J/T in CQM. In this model 
the hard Pomeron interacts with the constituent quark as it is shown in Fig. l-b. The vertex of the hard 
Pomeron for this interaction is 


V P ^3q = j\{d\* N{m) e qxi T fjf (rj), 


( 2 . 20 ) 


where 'k(rj) is the plane wave for the quark, and 'kjv is the wave function of the nucleon. This wave 
function actually depends on two variables. For example rq ^ = n — rq and rqq .2 = rq — \{r\ + V 2 )- The 
elastic differential cross section is proportional to 


\V p ^ 3q (t)\ 2 = /d 2 n,2dV3,l,2|^K2,r3,l,2)| 2 e— ^ 

1 + <? 2 y d 2 ri^d 2 r 3tl: 2r^ 2 \^ N ( r i,2,r3,i,2)\ 2 + 0(<? 4 ) = 1 + q 2 < |r 2 2 | > + 0(9 4 )- 


( 2 . 21 ) 


The value of < \r\ 2 | > can be evaluated from the q = y/—t behaviour of the elastic diffraction (see Fig. |j). 
This diagram leads to Ah ( t ) of Eq. ([2.16|) given in the form 


A H {t) = /d 2 n > 2dV3,i > 2|^(n, 2 ,r 3> i,2)| 2 ^3-3,i,2 


1 + - < |r 2 2 | > + 0(q 4 ). (2.22) 

Therefore, Bgg = § < |r 2 2 | > • Adding the combinatorial factors, we see that the t-dependence of the 
differential cross section for diffractive J/\k production is given by 


da(j* +p-> J/’k + 3q) 
dt 


oc 6 exp (~B G g \t\j + 3 exp (- B l £ G |i|) , (2.23) 


where the second term originates from the elastic rescattering off one quark. Using the same procedure to 
extract the slope B l GG from B in in Fig. [I], we estimate B GG = 1 GeV~ 2 . 

Eq. Q2.18|) has, thus, the form 


da H 

d?b 


1 


it 2 B G g 


exp 


b 2 \ 4 

+ 4 


1 


2 B gg 

1 


n 2 B l GG 


vr ( B GG + B% g ) 
b 2 

2 B%g 


exp 


b 2 


Bgg + B ; 


m 

GG 


+ 


exp 


(2.24) 


























2.4 Modelling the impact parameter dependence of the soft scattering amplitude 


In the eikonal model the opacity fl corresponds to a single soft Pomeron exchange. The t-dependence of 
this contribution is 


V#(i) 


So 


A - a' p |t| 


(2.25) 


where r(t) is the proton-soft Pomeron vertex with the normalization V(t = 0) = 1. 

We will use two simple models for the t-dependence of the vertex Vjy(t). 

1) An exponential parametrization valid in the forward t-cone. 

Bn 


V N (t) = exp - 


3 0 ,el 


(2.26) 


where Bq^i is the exponential slope of the elastic differential cross section due to a soft Pomeron exchange 
at s = so- Although there is no theoretical justification for it, the parametrization suggested in Eq. ( |2.26 ), 
is in accord with all experimental observables, that are sensitive to the small t region. In the following it 
will be denoted GP. 


2) A parametrization more suitable to cover a wider t range, including the diffractive dips observed in -^ L 

is 

Vjv(i) = G N (t). 

G]y(t ) is the proton electromagnetic form factor (see for example Ref. 

G N (t) = 1 


eterization for this form factor 


(1 + t/m 2 ) 2 ’ 

with m? = 0.72 GeV 2 . The above parametrization is used in CQM. 
The profile functions corresponding to GP and PP have the form 


(2.27) 

We will use the dipole param- 

(2.28) 


Scp(s,b ) = 


1 


2vr B, 


el 


exp 


b 2 


2 B, 


el 


(2.29) 


with B e i = B 0 e i + 2 a’p ln(s/s 0 ). 

Spp(s, b) = J qdqJ 0 (qb)G 2 N {t) exp (- a’ P ln(s/s 0 ) q 2 ) ■ (2.30) 

2.5 Explicit analytic formulae for the Gaussian parametrization 

Using the parameterization of Eq. (1.1), Eq. 0), Eq. ( |1.3| ) and Eq. ( |2.29| )) for the soft profile function, 
we can obtain simple analytic expressions for the main observables of the soft interaction || . 


o' tot = 2 J d 2 b{ 1 - e °0’ ft )/ 2 ) = 4vrR eZ (s) 
(Tin = J d 2 b( 1 

&el = ®tot &in ~ 27rB e i(s) 


,n(dd )+C -Ei(-dd : 

- a(s ’ b) ) = 2VT B e i(s) \ln(v(s)) + C- Ei(-u(s ))] , 


Zn(^) + C + Ei(-v(s)) - 2 Ei(-^) 
(Tel/vtot = h \ln(^r-) + C + Ei(—i/(s)) - 2Ei(-4fj\ / + C- Ei {-^) 


(2.31) 

(2.32) 

(2.33) 

(2.34) 
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where 


v(s) = Q(s,b = 0 ), 


(2.35) 


Ei(x) = \dt and C = 0.5773. From Eq. ( 2.34 ) we note that the ratio of (Tel/vtot depends only on the 
value of u(s) and, therefore, provides a possibility to find the value of v directly from the experimental 


The SP, given by Eq. CD, can also be calculated analytically if we assume Eq. ( |2.18[ ) for dan(b) / d 2 b. 
We obtain 

12 \ _ a(s)^(a(s),v(s)) 


sn = 


(2.36) 


The incomplete gamma function is 7 ( 0 , 2 ;) = fj z a ~ l e~ z dz and a(s ) = For the profile function of 

Eq. ( 2.30|) we cannot provide a set of analytic formulae. We have used the above formulae as a check of 
our methods and numerical computations. 


2.6 Constituent quark model 


The characteristics of the soft processes can be alternatively calculated also with CQM|3(|. In CQM we 
assume that all hadrons consist of constituent quarks (two for a meson and three for a baryon), and all 
scattering processes go through the interactions of the constituent quarks. In this model the state of three 
free quarks diagonalizes the interaction matrix, the proton state, and the diffraction dissociation state, 
which can be viewed as the expansion of the three quark plane wave function 

3 

Yl V q {ri ) = a 'kp + /3T d , (2.37) 

2—1 


with a normalization a 2 + (3 2 = 1. The constituent quark model can be viewed, thus, as a particular 
realization of a two channel model which we will consider below on a more general basis. 

For a constituent quark-quark scattering amplitude we use the eikonal formulae Eq. ( |Tt| ) and Eq. ([0|). 
Accordingly, the amplitude for the scattering of a quark i\ from the first hadron with a quark in from the 
second hadron is 


Ul,*2V 


(s, b) = i ( 1 — exp I — 




The opacity , 2 is given by 


n CQ M = 
* 1,*2 


(2.38) 


(2.39) 


co 


7T (8 B l g G + 4 a' p ln(s/s 0 )) 


s 

so 


exp 




8 Bq G + 4 ap ln(s/so) 


where b il}i2 = r iltkl ,h = r h - |(r fcl +n 1 ). 

In this paper we simplify Eq. ( |2.39| ) assuming that b ~ It is certainly correct at 

very high energies, since b 2 oc Aa' P ln(s/so) ^ ^ . More advanced calculations in the framework of 

CQM have been presented in Ref. [[3C|1 . 
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Having defined the opacities, we obtain the following formulae for the soft cross sections 


Vtot 

d&ei/dt 


d&diff/dt + da e i/dt 


9x2 J d 2 b (1 - exp(-in c « M (s,&))) , 

9 2 x G%(t)\f c ^ M (s,t)\ 2 , 

9 2 x (|exp [B GG - B% a ] |t|) + 1 ) 2 \f CQM (s,t )\ 2 . 


(2.40) 

(2.41) 

(2.42) 


fCQ M (s,t) is the quark-quark scattering amplitude in the t-representation, see Eq. ( ^H) and Eq. 
for the relations between t and b representations. Gjy(t) is the electro-magnetic proton form factor (see 
Eq. (2.25)). The numerical combinatorial factors give the number of interacting quark-quark (antiquark) 
pairs. Note that in the CQM diffraction is confined to small masses] 

In CQM, with d 2 aH/d 2 b given by Eq. ( 2.24 ), we obtain 


< \S 2 \ > = 


f d 2 b ( da H /d 2 b ) e 9 ' n H.i 2 <>> b ) 
f d 2 b (d<JH/d 2 b) 


(2.43) 


where is determined by Eq. ( 2.39| ). Eq. ( 2.43| ) determines the probability that none of the possible 

nine quark-quark pairs can rescatter inelastically leading to the production of additional hadrons. 

Comparing the calculations using this model with the calculations using the eikonal approach, we are 
able to assess how sensitive the value of SP is to a more elaborate model. 


2.7 Two channel model 

A major deficiency of the single channel eikonal model is its failure to reproduce the observed mild energy 
dependence of a sc i in the ISR and above||. In a single channel model, cr sc i is assumed to be small enough 
so that is a small parameter. Hence, the eikonal rescatterings are elastic. Accordingly, <j s d is not 
included in the fitted data base fixing the single channel model parameters. Regardless, the SP of a s d can 
be calculated resulting with disappointing output. The model’s inability to reproduce the diffractive cross 
section energy dependence reflects an over estimation of SP at higher energies (see Ref. Q for details). 

To overcome this difficulty we have developed a more elaborate multi channel eikonal model0 in which 
both elastic and diffractive soft rescatterings are included resulting in a decrease of the calculated high 
energy SP. In the two channel version of the model double diffraction is assumed to be small enough so 
as to be neglected in the eikonal rescatterings. This assumption is supported by the data available in the 
ISR-Tevatron range. Note that the addition of a competing diffractive channel actually reduces the elastic 
screening. In our model, the added diffractive screening adds more screening than the loss it has induced. 
The final results in an over all larger screening which reproduces the SD data well, considering the scatter 
of the experimental points |24fl . 

In a two channel model, the two hadron states of a proton and a diffractive state are considered 
simultaneously. 


= afi + (5^2 , (2.44) 

= -/3Ti + a^ 2 • (2.45) 
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N 



Figure 5: Production of di-jet with two large rapidity gaps in the two channel eikonal model for the soft interaction. 


p and D denote the proton (Eq. (|2.44|) ) and diffractive (Eq. (|2.45| )) states. The wave functions functions 
Ti and 4^2 are orthogonal. They diagonalize the interaction matrix at high energies and, therefore, 


— Ai,k $i,j $k,l ■ 


As in the single channel formulation 


d 7. A' ^ I 1 


iO 

_ p 2 ^ “i,k 


For the opacities 12^(6) we use a GP parameterization 

- b2 

^ i,k ( b ) = ^,fc e ^ = 

where <rg >fc and 




* R lk( 8 ) V«o 


. A b 2 

S \ 

— e i ’ k , 


R i,k( s ) ~ % R i ,o + 2^,0 + 4 ap ln(s/so). 


(2.46) 


(2.47) 


(2.48) 


(2.49) 


More details regarding the parametrization can be found in Ref. Jc]. An important fitted parameter which 
is essential for the present study is /? = 0.464. 

For the calculation of exclusive di-jet central production SP we need to generalize Eq. (|L^) for a two 
channel scenario. This generalization is illustrated in Fig. |5], leading to a generalization of Eq. (2.8) 


A (N + N -*• N + LRG +JJ+ LRG + N) = 


(2.50) 
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O'h.ard ^ ' 
i,k 


d 2 b\ d 2 b 2 e ipu + * 3?2t '^ 2 An(i —> p; &i) 7l//(fc —»■ p; 62 ) e' 


»»,fc(( g l+ g 2) 2 


This equation can be rewritten in an explicit way using Eq. ( 2.44| ) (see Refs.[[7], 241 for details). 

A (AT + AT AT + Li?G + JJ + LRG + JV) = (2.51) 

= a hard J d 2 b 1 d 2 b 2 e ipit ^ 1 + ip2t ^a 2 e ^ “ 4 ff(l -> p; 61 ) A H (l p; b 2 ) 

-ni, 2 ((bi -1-62) 2 ) 

+ a/3e 2 {A h (1 p\bi) A h {2 -^p; 6 2 ) + A H {2 -> p; 61 ) R#(l -> P] h 2 )} 

~ n 2,2((*>l +^2) 2 ) \ 

+ (3 2 e 2 A H (2 ^ p-,bi) A H (2 ^ p-b 2 ) . 


Eq. ( 2.51 ) can be written in a more convenient form which enables us to use the experimental data 
corresponding to Fig. f|. We introduce two amplitudes 

V p ^ p ( b 2 


A H (p —> p; b) = 


2 ttB gg 


exp 


Ah(p —> D\b) = exp 

Z7ri3 GG 


2 B gg 
b 2 

"2 b% g 


(2.52) 

(2.53) 


where the input parameters Bqq and B l GG have been introduced in Section II.C. The input assumption of 
the two channel model[j7j] is that the double diffractive production is small enough to be neglected. Using 
this assumption, together with Eq. (2.52|), Eq. (2.53|), Eq. (2.44) and Eq. (2.45|), enables us to rewrite 


Eq. (R5l|) 


A(N + N -> N + LRG +JJ+ LRG + N) = 

-fif(6i + b 2 f 


(2.54) 


/ 


= CLhard / d b\ d b 2 e 


exp 


1 -2 P 2 s d [(bi + b 2 ) 2 )) A H (p —> p;bi) A H {p -»■ p;b 2 ) 

- 2a/3 s D [(bi + & 2 ) 2 ) {A H (pp-,b{) A H (p-> D;b 2 ) + A H {p ->■ D\b{) A H (p -> p; & 2 
In our notation 

, h , , ( a n(b)\ 

s D {b) = 1 - exp I--— 1 , 

with A Q(b) = n h2 — Oip. 

For the opacities 12 and AQ we use the forms and fitted parameter values of Ref. @. 


(2.55) 


a 0,p j 

(s) 

1 exp | 

( *> 2 ^ 

7tR 2( s ) 

l S oJ 

l WAJ 

a ° ,D | 

(s) 

A 

1 exp | 

( 62 

TTRjj(s) 


\ Rd( s ) ) 


(2.56) 

(2.57) 

(2.58) 
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where 


Note that Rq D = \Rq p - 


RporD — Rq, porD + 4 Up lll(s/so). 


To obtain the SP we need to substitute the following two elements in Eq. (1.5) 

d&H r^H 


d 2 b 


TS(b), 


and 


—Cl dcrH 
3 d 2 b 


e- Q ((l-2/3 2 sp) 2 T e f(6) 


- Ava[3(l-2 0 2 s D )s D T?(b) + 8 v 2 a 2 (3 2 s 2 D { if (b) + if (5)}) 
v = V p ^ D /V p ^p, and 

1 


rpH _ 

1 el ~ 


rri JJ _ 

T i = ZTRh ex P 


2 IT B gg 

1 


exp 1 -db) 

b 2 


The corresponding input slopes are 


B? 

? H 


irBj 1 " V B?) ' 
1 + 3k 


1 + K 


Bgg , 


b 2 = (1 + «) Bgg , 


B? = 


4 K 


1 + K 


Bgg , 


where n = B 1 £ g /B G g- 


3. Results and Comparisons 


(2.59) 


(2.60) 

(2.61) 


(2.62) 

(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 


In this section we present the results of the four models we have considered and compare them with the 
relevant experimental data. We then proceed to present and assess our SP predictions aiming to define a 
value and a range for the SP of a central GJJG exclusive dijet production at the LHC. 


3.1 Adjusted parameters 

Following is a summary of the phenomenological adjusted parameters for the soft input of the four models 
we have considered. 


The intercept of the soft Pomeron trajectory at t = 0 

A, 

(3.1) 

The slope of the soft Pomeron trajectory at t = 0 

a P , 

(3.2) 

The initial energy squared 

so, 

(3.3) 

The strength of the Pomeron interaction at s = so 

°0 j 

(3.4) 

The slope of the vertex Vi\r(t) 

Bo,el , 

(3.5) 

The slope of the vertices in the two channel model 

R-o,p an( 4 Ro,d i 

(3.6) 

The fraction of Tp and 4*2 

/5• 

(3.7) 
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We determine these parameters by fitting the value and energy dependence of the soft data base observables. 
As previously mentioned, the input assumption of a single channel eikonal model is that « 1. 

Accordingly, its data base consists of atot, &ei and B e i, whereas a s d is a prediction. CQM parameters are 
adjusted from the same data base. Since this model applies only to small mass diffraction^], it cannot 
predict a s d■ The data base for the two channel model includes the above and in addition the a s d data 
points. We calculate the corresponding SP as a prediction derived after fixing these parameters, provided 
we can specify the opacity (opacities) of the screened hard process. 

The best fit adjusted parameters for the models considered are presented in Table I. Note that even 
though the fitted ctq values obtained for the single channel model GP and PP profiles are identical, the 
corresponding v values are different reflecting different b-distributions. 

3.2 Reproduction of the soft scattering observables 

A quality reproduction of atot, &ei and B e i is an obvious prerequisite for a soft scattering model with which 
we may calculate the SP of interest. As we noted, a s d is not included in the soft data base of the single 
channel eikonal model, but it is a prediction of the model IQ]. CQM is not suitable to calculate a s d ■ In the 
two channel model, a s d is included in the fitted data base. In the following we discuss the details of the 
soft scattering output of the models we have considered. 

3.2.1 One channel model 

For the purpose of our present investigation we have considered a toy soft DL like||22|] Pomeron exchange 
model, neglecting the secondary Regge exchanges. This is a strictly high energy model for which we take, 
never the less, a relatively low sq = 450 GeV 2 . This choice should be compared with standard Regge 
pole parametrizations in which the Regge contribution diminishes at much higher energies. Consequently, 
this model can not be continued to lower energies. The best adjusted Pomeron parameters are given in 
Table I. Regardless of its simplicity, this model provides a very reasonable reproduction of its data base 
in the ISR-Tevatron enery range. An interesting observation is that the two profiles result in remarkably 
close outputs. This is so even though the v values corresponding to the two profiles are quite different. 
i/ gp ( 450) = 1.88, while vpp( 450) = 2.31. We conclude that a good reproduction of the data base requires 
a delicate compensation between v and the effective radius of the two profiles respectively. 

Checking the b-distribution output of a given eikonal model, we encounter a fundamental feature, 
which is very transparent in the single channel version with a GP profile input for both the elastic and 
diffractive amplitudes. The elastic amplitude output, which respects unitarity, maintains an approximate 
Gaussian b-distribution peaking at b = 0. On the other hand, the diffractive amplitude output is non 
Gaussian peaking at a non zero b. The sum of the two amplitudes respects the Pumplin bound. This 
phenomena is easily understood once we note that the eikonalization of a diffractive amplitude amounts 
to the transition Mdiff(s,b ) —> Mdiff(s,b)e~ n ( s ’ b \ Since Q is central its suppressing effect is maximal at 
small b, and diminishes at high b (For details see Ref. II and the preprint first version of this paper |]3l| ). 

As noted, the problematic asset of the single channel model is its inability to reproduce the energy 
dependence of a s d, which implies an over estimation of the corresponding SP becoming worse with increasing 
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Figure 6: CQM LHC prediction for the t-dependence of the low mass diffractive differential cross section. 

energy. As a result, we consider any SP calculated within the single channel to represent, at best, an upper 
limit for the SP. 

3.2.2 CQM 

The CQM does not provide a good reproduction of the soft data base at low and medium energies. It is 
applicable for the Tevatron energies and above. Note that in this model the PP profile is naturally related 
the electro-magnetic form factor of the proton. CQM is limited to small mass diffraction only, since the 
triple Pomeron vertex is not included in its formulation. Consequently, a s d is not included in its soft data 
base. Since the model takes into account the rescatterings due to low mass exitation of the interacting 
hadrons it may be classified as an example of a simple two channel model. 

The advantages of the model is it simplicity and the realistic b-dependence. The model does enable us 
a calculation of low mass drTd j/ f /ctdif f presented in Fig. |b| at LHC energy. The positions of the calculated 
dips are a consequence of the input PP profile and should be considered reasonably reliable since the model, 
with this input, correctly reproduces the t dependence of the Tevatron elastic cross sections. 

3.2.3 Two channel model 

Our reproduction of a tot, &eh B e i and a s d in a two channel eikonal model have been published a while agoj/J]. 
Note that the soft data base includes the single diffraction data. Our d x o j = 1.5, which is seemingly high, 
reflects the poor quality of the a s d points. Note that the b-space peripheral behaviour of the diffractive 
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channels amplitudes is maintained in the two channel scenario. Our LHC predictions for crtot = 103.8 mb, 
u e i = 24.5 mb, B e i = 20.5 GeV~ 2 and a S( j, = 12.0 mb. We consider these predictions, and the consequent SP 
calculated values, reliable as our input diffractive rescatterings are based on an effective parametrization 
which includes the high mass exitations. The above numbers are significant for central Higgs diffractive 
production, where there is a problematic background of single diffraction which mimics the sought after 
Higgs signal. 


3.3 Survival probabilities 


Using Eq. (|1.5|) we calculate the energy dependence of the SP for our process in the single channel eikonal 
model. The result is shown in Fig. [?], which compares the exclusive central di-jet SP as calculated in the 
four models we have considered. For a single channel model we see that the profile function for the soft 
interaction does not considerably affect the value of the SP. The GP and PP profiles produce a small SP 
difference at the ISR energies, but at higher energies, Tevatron and above, the results are essentially the 
same. 


In Section II.C we have explored the two component structure of ypf-, (Eq.Eq. (2.24)). Fig. ^examines 
the impact of the inelastic component in the calculation of SP for our process in a single channel model 
with a PP profile. The results show a relatively small difference (within 10% accuracy) between a single 
elastic component and the sum of elastic and non elastic components (see II.C for details). An increase of 
the fraction of the inelastic production does not change the value of the SP significantly. 


Fig.0 details the SP predictions of the four models we have examined. We aim, on the basis of these 
predictions, to suggest upper and lower bounds for the SP corresponding to exclusive central diffractive 
di-jet production at the LHC. Examining Fig. |we observe that SP calculated in the two channel model is 
consistently lower than the corresponding single channel model and CQM predicted SP values. This may 
serve as a guide in our attempt to suggest a margin of error in the determination of SP at the LHC. 


When assessing to higher SP bound, we note that the CQM SP is much lower than the single channel 
values at W < 1000 GeV, with a difference that gets smaller with increased energy. The two predictions 
are approximately the same, 6.0±0.1%, at the LHC energy and they cross just above the LHC. This 
suggests a 5-6% as the upper bound for the calculated SP at the LHC. We consider the SP estimates with 
models neglecting the diffractive channel rescatterings to over-estimate the calculated SP output. Prudency 
suggests, thus, an upper SP bound of 4- 5%, which is moderately smaller than the predictions of the single 
channel models and CQM. 


The two channel eikonal prediction for exclusive central diffractive di-jet production at the LHC is 
2.7%, compared with 3.6% for the corresponding inclusive central di-jet production^]. The two channel 
input is v = \/3, Bgg = 3.6 GeV~ 2 and B GG = 1.0 GeV~ 2 . Our two channel predicted SP at the LHC is 
almost identical to the KKMR[32] value. This is very supportive of a SP = 2.5-3.0% at the LHC. 


It is instructive to compare our model with the KKMR model[^]. The two models are defined as ’’two 
channels”, but are, actually, rather different. 

1) Our two channel eikonal definition, for either soft or hard diffraction, consistently refers to two possible 
modes of soft rescattering, i.e. elastic and diffractive. Accordingly, both elastic and diffractive states in 
a p — p scattering are presented as a linear combination of our two orthogonal base wave functions VR i 
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Figure 7: Energy dependence of the survival probability for exclusive central di-jet production as calculated with the 
four models we have considered. 


and 'J / 2 - In our two channel model we neglect the double diffraction channel, which is exceedingly small 
in the ISR-Tevatron energy range. The screening opacities of our input eikonal matrix have different b=0 
normalizations, i.e. different v values, and different b-dependences which reflect our different input for the 
elastic and diffractive forward differential slopes. 

2) In the KKMR model for soft interactions the eikonal matrix is defined in a similar way to ours. Note, 
though, that their diffractive eikonal components are restricted to low diffractive mass. The screening 
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Figure 8: Energy dependence of the survival prohahility for di-jet production in the eikonal single channel model with 
a PP profile. The upper curve corresponds to the final state production of two nucleons. The lower curves 
take into account the inelastic production of excited states assuming either Eg. ( 2.2^ ), or that at t=0 the 
elastic and inelastic production have the same amplitude. The two last cases can hardly be distinguished. 


opacities have the same b-dependence, identical to the b-dependence of the single channel elastic opacity, 
having a different b=0 normalizations. Unlike our input, KKMR do not neglect the double diffractive 
channel. 

3) It is no surprise that KKMR obtain a very high LHC prediction for a^d = 9.5 mb, which is comparable to 
their predicted a s d. KKMR predict, thus, an inelastic diffractive cross section a s d + a^d = 18.9 — 24.9 mb, 
relative to a atot = 99-1 — 104.5 mb. Our predction is a s d = 12.0 mb, which may be corrected by a small 
value add, , relative to a atot = 103.8 mb. Note that we were not successful in reproducing the KKMR fit for 
a e i and B e i. KKMR have not published a detailed prediction for a s d- Their quoted LHC value for single 
diffraction is comparable to ours. 

4) In the KKMR model the hard process has two dynamical components treated as two eikonal chennels. 
Two sets, with two components each, were studied and have been shown to produce very similar results. 
In our two channel model the hard diffractive proces has two components associated with the possibility of 
either an elastic or inelastic diffractive initial rescattering preceeding the hard process. The SP is calculated 
accordingly. 

5) In our model we make a weak assumption in which the ratio between the elastic and diffractive couplings 
is the same for soft and hard Pomerons. Consequently, we are able to associate the hard amplitude with 
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HERA DIS experimental data resulting in eikonal opacities which are different in their b=0 normalization 
as well as their b-dependence. KKMR make a much stronger assumption in which all opacities, soft and 
hard, have the same b-dependence and a normalization which is determind by the above coupling ratio 
assumption. 

6 ) Regardless of the above the actual KKMR SP results are remarkably close to ours. This result deserves 
a clarification in the future. 


4. Transverse momentum dependence of the cross section 


In the following investigation we have calculated the dependence of the output di-jet cross section on pi d 
and p 2 ,t- This calculation provides the differential information on the t dependence of the cross section 
and SP of the process. The importance of this differential calculation is that it offers a refined method 
to discriminate between different models and/or parametrizations leading to compatible values of SP. We 
define a damping factor 


Dsp = 


_ (id _ 

d 2 pit d 2 p2t dyi < 4/2 


Id 2 Pitdip 2 t , i[ 2 pit d/ p l 4yidy2 


In a single channel eikonal model, 


Die = 


f d 2 bid 2 b 2 e i ™^ +i P*^A H (b 1 )A H {b 2 ) A s (( 6 i + b 2 y 
f d 2 bid 2 b 2 \A s ((bi + 6 2 ) 2 ) A H {bi) A H {b 2 )\ 2 


2 \ 2 


(4.1) 


(4.2) 


From its definition the value of the SP is dependent on the b-distribution of the soft amplitude profile. 
Should we have an experimental information on the dependence of Dgp on pi t and p 2d , we hope to be 
able to invert the above procedure and obtain information on the form of the b-profile from the differential 
properties of the SP. 


For the generalization to a two channel model we insert Eqs. 2.54 - 2.65 into the expression for the 
two channel SP. Using the amplitude, given by Eq. ( [2.54 ), we can carry out the integration over d 2 bi 
analytically, obtaining a much simpler expression. 


D 2 c = 


w 


f d?pi )t d 2 p 2) t\A\ 2 ’ 

where A = Ai + A 2 + A 3 . Defining b = b\ + b 2 and p\ 2 = p\ d — p 2 j, we obtain 
1 


Ai = 


4 PBl 


J d 2 be 2 ) s( 6 )exp^— 


b 2 1 9 \ (b 

J Jo [jIPP + P 2 .t 


el 


Ao = — 


2a (3 


, 2 , 


4vr 2 B P , B 2 


d~b e 


■ exp 


e l a inel • 
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sp(b) 


2 (Bi ne i + B e i 

A 3 = A 2 (pij ^ p 2jt ). 


{6' + Bi ne i B e ip \ 2 } ) To 
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inel 
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el 


\DelPl,t T Bi ne i P 2 ,t 
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(4.4) 

(4.5) 
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Fig. U and Fig. 0 show the single channel di-jet factor D\c with GP and PP b-profiles as a function 
of the transverse momenta of the produced protons in the final state. We predict a typical minima at 
pi t ~ 1.75 — 2 GeV. This behaviour is a manifestation of the wave nature of our diffractive scattering 
process. It depends on the scale of the soft profile b-distribution. From Fig. || and Fig. 0 one can see 
that the minima pt position depends on both the energy of p\.t and P 2 j, and the angle 9 between them. 
As stated, we hope that a future measurement of this dependence will provide valuable information which 
will add to our knowledge of the soft interaction amplitude in QCD. To get some initial information on 
the angular distribution, we have chosen (for convenience) 6 = 0 , 7r/2,and7r. 


In Fig. [Tl] and Fig. 12 we compare the transverse momentum distribution of the di-jet production 
differential cross section with the corresponding elastic cross section. These figures illustrate that the di-jet 
cross section has quite a different structure of minima than the elastic cross section, where the positions 
of these minima move to smaller values of pt at higher energies. Comparing with the experimental data 
(for data see Ref. 2£] and references therein), we deduce that this data is compatible with the PP profile 
describing the structure of the minima and maxima in proton-proton collisions in the energy range of 
W = \/s = 20 — 65 GeV. Our parametrizations for the profile function give a range of predictions for 
the di-jet production at the LHC. The difference between the GP and PP parametrizations, as far as pt 
dependence of the di-jet cross section is concerned, is much smaller than for the elastic cross section. 


A most interesting result is the striking difference between the t-dependence of the damping factors, 
defined in Eq. obtained in the single and two channel eikonal models. This difference is clearly seen 

in Fig. 13 in which GP profiles were used. The single channel model leads to a dip and to a slope at pi,i=0 
which is equal to Bqp without a suppression due to SP. In the two channel eikonal model we do not have 
any dip at p± t t < 2 GeV and the slope turns out to be much smaller than the slope obtained from the 
single channel eikonal model. 

This result can be anticipated from the general formula of Eq. ( |2.51 ) and Eq. ( [2.53 ) which show that 
the damping factor is proportional to a contribution of the order of /3 while the elastic cross section is 
proportional to /3 2 (see Ref.j7|). 


5. Conclusions 

In this paper we have calculated the large rapidity gap survival probability < \S\ 2 >, for exclusive central 
dijet production at the LHC. Our assessments are of particular interest as the SP we have calculated 
correspond also to a central exclusive diffractive production of Higgs. The unique feature of our treatment 
is that we explicitly included the impact parameter dependence of the hard amplitude in our evaluation 
of SP. This b-dependence appears in the proton-hard Pomeron vertex, whose structure was deduced from 
information obtained from the DIS production of J/il>. 

The calculation of the LRG SP requires knowledge of both the hard and soft components of the 
hadron-hadron interactions. As there is no reliable theory for the soft component, we employed four 
phenomenological eikonal models to describe the soft processes. These models are: 

1) A single channel model with an exponential dependence in t for the proton-soft Pomeron vertex. 

2) A single channel model with a power-like behaviour for the proton-soft Pomeron vertex. 
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3) The CQM, which is a non Regge two component soft model, where we have assumed that this vertex is 
given by the proton electromagnetic form factor approximated by a dipole. 

For cases 2) and 3), we have to integrate over b numerically, while for case 1), we have an analytical 
expression. 

4) To these we add the two channel model whose main advantage is that it is the only model considered 
that properly reproduces the available diffractive data. 

All four models considered were utilized to calculate the SP of an exclusive central hard di-jet produc¬ 
tion. The two channel model predicts for this SP a value of 0.027, which is almost identical to the KKMR 
prediction. These should be compared with a 0.036 value for the equivalent inclusive calculation |24(|. The 
other three models have compatible prediction of 0.06. We consider this value to be some what over es¬ 
timated as these models neglect the diffractive rescattering which increases the eikonal overall screening 
and, thus, reduces the SP. We prudently conclude that the LHC SP we have studied lies between 0.03 - 
0.05. 

The single channel cross-section for the dijet production as a function of pt, has a typical structure of 
minima at certain values of pt, which is dependent both on the angle between p t \ and pt 2 of the two jets, 
and on their energy. Positions of the minima are model dependent. 

The most interesting result of this paper is the fact that the two channel model predicts quite different 
dependence of the damping factor versus the transverse momenta. The central prediction is that there 
are no dips for pt < 2 GeV. Also, the slope at pt = 0 turns out to be much smaller than for the one 
channel eikonal model. Both features stem from the general properties of the two channel model and we 
firmly believe that this prediction should assist to select the appropriate models for the description of SP 
directly from the experimental data. Since the general behaviour of the damping factor depends weakly 
on energy we suggest assessing the damping factor at the Tevatron energies. We would like to stress that 
the calculation of the damping factor show such simplicity and elegance that we believe they have a deeper 
meaning than being just a consequence of a particular model. 
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Parameters 

A 

a' P { 0) 

so 

co 

Slope 

Gaussian parameterization 
(GP, Eq. (2.29)) 

0.09 

0.25 GeV -2 

450 GeV 2 

47.2 mb 

B 0 ,ei= 10.24 GeV 2 

Power -like parameterization 
(PP, Eq. (2.30)) 

0.09 

0.25 GeV” 2 

450 GeV 2 

47.2 mb 

m 2 = 0.72 GeV" 2 

Constituent Quark Model 
(CQM, Eq. (2.39)-Eq. (2.42)) 

0.08 

0.28 GeV -2 

250 GeV 2 

4.13 mb 

B% g = 0.5 GeV -2 

Two channel model 

2ChM 0 = 0.464 
(Eq. (2.44) - Eq. (2.48)) 

0.126 

0.2 GeV~ 2 

1 GeV 2 

5.07 mb (<Tp) 
56.5 mb (cr%) 

16.34 GeV -2 (i$ iP ) 
8.17 GeV" 2 (R\ d ) 


Table 1: The best fit adjusted values of the Pomeron parameters in the models considered. 



















Figure 9: Transverse momentum dependence of the cross section factor D in a single channel model with a Gaussian 
b-profile (GP) for various 9 values. We also show an exp(—3.6|£|) dependence corresponding to to the 
hard slope Bqg- 
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Fig- [lOf-e 



Fig. lC-c 



Fig. [l^-b 


Figure 10: Transverse momentum dependence of the cross section factor D in a single channel model with a power 
like h-profile (PP) for various 9 values. We also show an exp(—3.6|t|) dependence corresponding to to 
the hard slope Bqg- 
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Fig. [TTJ-a 



p,, [GeV] 


Fig. IK-b 


Figure 11: Transverse momentum dependence of the cross section factor D in a single channel model with a Gaus¬ 
sian b-profile (GP) compared with the corresponding elastic cross section. 
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Fig. 12 


Fig. 12-b 


Figure 12: Transverse momentum dependence of the cross section factor D in a single channel model with a power 
like b-profile (PP) compared with the corresponding elastic cross section. 
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Fig. 13-a 



Fig. 13-c 



Fig. [l3|-b 


Figure 13: Transverse momentum dependence of the cross section factor D in a two channel model with Gaussian 
b-profiles (GP) for various 0 values, compared with the corresponding single channel elastic cross section 
with 6 = 7t/2. We also show an exp(—3.6|t|) dependence corresponding to to the two channel input hard 
slope Bgg- 
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